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IX 


1.  Objective 


The  objective  of  this  study  was  to  develop  a  user  friendly  model  to  investigate  heavily 
loaded  lubricated  circular  and  noncircular  (elliptic)  contacts.  In  order  to  achieve  the  objectives, 
a  numerical  model  for  the  Newtonian  thermal  elastohydrodynamic  lubrication  (EHD)  of 
rolling/sliding  point  contacts  was  developed.  The  model  simultaneously  solves  the  two- 
dimensional  Reynolds,  Elasticity  and  the  three-dimensional  energy  equations  to  obtain  the 
pressure,  film  thickness  and  temperature  distribution  within  the  lubricant  film.  The  numerical 
model  was  extended  with  an  interactive  computer  graphics  to  aid  in  the  design  and  analysis  of 
elastohydrodynamic  lubrication  of  rolling/sliding  tribo-contacts.  This  report  provides  a 
description  of  the  method  of  solution,  the  computer  aided  design  graphics  package  and  a  sample 
of  the  results. 
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2.  Introduction 


The  lubrication  of  heavily  loaded  high  speed  machine  components  such  as  bearings,  gears, 
etc.,  has  been  extensively  investigated  in  recent  years.  However,  most  of  the  analyses  have 
routinely  assumed  the  lubrication  process  to  be  isothermal.  Most  of  the  research  and  published 
works  have  dealt  with  isothermal  and  thermal  EHD  lubrication  of  line  contacts  or  isothermal 
EHD  of  point  contacts.  The  study  of  heavily  loaded  lubricated  line  contacts  with  the  thermal 
effects  included  started  with  Stemlicht  et  al.  in  1961.  Cheng  (1965)  investigated  the  thermal 
effects  in  EHD  lubrication  of  »»ne  contact,  however,  the  loads  considered  in  his  analysis  were 
extremely  low  which  resulted  in  insignificant  temperature  rise.  Murch  and  Wilson  (1975) 
illustrated  that  temperature  wi.'  significantly  affect  the  film  thickness  particularly  at  high  speeds. 
Sadeghi  and  Dow  (1987)  used  experimental  pressure  and  surface  temperature  measurements  to 
investigate  the  temperature  rise  within  the  lubricant  film  operating  under  EHD  lubrication.  They 
illustrated  that  the  temperature  rise  is  significant  and  cannot  be  neglected.  Recently,  Sadeghi 
and  Sui  (1990)  developed  a  complete  numerical  solution  to  the  EHD  lubrication  of  line  contacts. 

The  EHD  lubrication  of  point  contacts  has  been  investigated  less  rigorously.  This  is 
mainly  due  to  the  extreme  complexity  involved  in  the  solution  of  the  two-dimensional  Reynolds 
and  elasticity  equations.  Many  investigators  have  studied  the  isothermal  EHD  lubrication  of 
point  contacts,  for  example.  Ranger,  Ettles  and  Cameron  (1975),  Hamrock  and  Dowson  (1977), 
and  Evans  and  Snidle  (1982).  However,  most  of  these  analyses  are  extremely  computer  time 
intensive.  Recently,  Lubrecht  et  al.  (1987)  used  multigrid  technique  to  solve  EHD  lubrication  of 
point  contacts.  This  method  provided  a  fast  converging  scheme,  however,  they  neglected  the 
thermal  effects  within  the  lubricant  film.  Briiggemann  and  Kollmann  (1982)  are  among  the  first 
to  include  the  thermal  effects  in  EHD  lubrication  of  point  contacts.  They  only  considered  the 
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conduction  across  the  film  and  viscous  dissipation.  They  chose  to  neglect  compression 
neating/cooling  and  convection  terms.  However,  it  has  been  shown  (Sadeghi  and  Sui,  1990) 
that  compression  heating/cooling  plays  a  major  role  in  EHD  lubrication  and  since  the  Peclet 
number  is  usually  high  the  convection  term  cannot  be  neglected.  They  also  assumed  a  pressure 
distribution  rather  than  solving  the  Reynolds  equation.  Zhu  and  Wen  (1984)  presented  a 
solution  for  thermal  EHD  lubrication  of  point  contacts.  However,  the  loads  considered  in  their 
analysis  were  extremely  low.  Blahey  and  Schneider  (1986)  also  developed  a  solution  for  the 
thermal  EHD  lubrication  of  elliptic  contacts  using  control  volume  approach.  However,  they  also 
neglected  the  convection  term  in  the  energy  equation.  Recently,  Faghri  (1984)  presented  a 
methodology  for  finite  difference  solutions  of  two-dimensional  convection-diffusion  problems  in 
irregular  domains,  using  a  non-orthogonal  coordinate  transformation. 

In  this  report,  the  graphics  package  for  the  EHD  analysis  and  the  technique  used  to  solve 
the  three-dimensional  energy  equation  for  thermal  EHD  lubrication  of  point  contact  under 
rolling/sliding  condition  is  described. 
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3.  Mathematical  Formulation 


Figure  1  illustrates  a  schematic  of  the  problem  under  consideration.  The  lower  body  (flat 
plate)  and  the  upper  body  (an  ellipsoid)  are  moving  at  the  surface  velocities  of  u,  and  u2, 
respectively.  Lubricant  is  pulled  into  the  contact  due  to  the  viscous  effects  of  the  lubricant  and 
the  movement  of  the  bodies.  The  load  causes  the  pressure  within  the  lubricant  film  to  elastically 
deform  the  bounding  surfaces. 


Reynolds  Equation 

The  dimensionless  Newtonian  thermal  Reynolds  equation  with  the  appropriate 
assumptions  derived  from  the  Navier-Stokes  equation  can  be  written  as  [Dowson,1962]: 


ax 


dY 


pH3|P  [Gl-C2]  +[X2-X1)  Jr[pHG3)+X1 


where 
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FiCD  =  1  -  p  T0  (T  -  1), 
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Ellipsoid 


Figure  1.  Lubricated  contact  of  an  ellipsoid  and  a  semi-infinite  plane. 


p  =  p  at  T=  1, 


(7) 


The  boundary  conditions  for  equation  (1)  are: 


P  =  0  on  the  periphery 
3P  3P 

P  =  ~  =  —  =o  on  the  cavitation  boundary 

oX  a  Y 


The  constant  load  condition  is: 

7  7  27tke 

J  J  P  dX  dY  =  — — 


(8) 

(9) 


(10) 


Energy  Equation 


The  Reynolds  equation  is  two-dimensional,  however,  the  energy  equation  is  three- 
dimensional  because  of  the  temperature  variation  in  the  Zt -direction  (Figure  2).  The  energy 
equation  in  dimensionless  form  is[Sherman,  1990]: 
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The  dimensionless  boundary  conditions  for  the  energy  equation  (11)  are  (Carslaw  and  Jaeger 
,1959  and  Zhu  and  Cheng,  1989): 
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(13) 
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Film  Thickness  Equation 


The  elastic  deformation  of  the  surfaces  in  dimensionless  form  is  given  by: 


X2  +  Y2/r  20 
H(X,Y)  =  H„  +  A  v-  '  +-5 


iZL  f  f  PdX'dY' 

n2  J  2  [(X-XO2  +  (Y-YO2]1 


jcPhR* 
where,  0!  =  £,-  — 


Viscosity-Pressure-Temperature  Relationship 


The  viscosity/pressurc  relationship  used  in  this  analysis  was  proposed  by  Roelands  (1963) 
and  further  developed  by  Houpert(1985).  The  Roelands  equation  in  dimensionless  form  is: 


q  =  exp 


[ln(rio)  +  9.67]  y  (l  +  5.1  x  KT9  PH  P]¥  - 


Density-Pressure-Temperature  Relationship 

The  density/pressure  relationship  used  by  Dowson  and  Higginson  (1966)  in  the 
dimensionless  form  is  employed  in  this  analysis. 


p=  1  + 


0.6  x  ior9PHP 
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Coordinate  Transformation 


The  energy  equation  (11)  requires  a  coordinate  transformation  due  to  the  curved  boundary 
before  control  volume  formulation  can  be  applied.  The  transformation  equations  used  in  this 
study  are: 


\-X,  S-Y,  and  C 
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The  unit  vectors  along  the  new  coordinates  are  ,  es,  and?;  where .  Figure  2a  illustrates 
the  unit  vector  n  normal  to  the  curved  plane  of  C  =  constant  at  a  point  Q: 


-  VC 

n=7vcr 


-  ^  -  7  Cy  +?z, 


where  a*  Vl  +  j52  +  y^,  P  =  and  y=  a£-|^~ 


The  unit  vector  n  is  perpendicular  to  the  line  of  constant  4  and  line  of  constant  6  passing 
through  the  point  Q,  therefore  the  unit  vectors  e\  and  e8  are, 


+  P  3Z, 


ViTjf 


_  ey+YCz, 
c«  =  — p=- 

Vl+71 


The  velocity  vector  in  the  new  coordinate  system  is 


U  =  U5e;  +  V8e»  +  W;e; 


and  the  velocity  components  are 


VTTfu 

V,  =  Vl  +  Y  V 


w?=w-jSu-yv 


The  partial  derivates  in  terms  of  transformed  coordinates  are 
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J_  ±  J_± 
ax  "  ac  -  ah  ac 


_d_  _  _3 _ Y_  _a_ 

ay  =  as  “  ah  ac 


a  i  a 
az,  =  ah  ac 


and  the  gradient  operator  V  is: 


V  = 


_a_  X  A 
ac  ~  ah  ac 


cx  + 


J_  jl  JL 
as  "  ah  ac 


eY  + 


i  _a_ 
ah  ac 


Control  Volume  Formulation  and  Discretization 
The  Reynolds  equation  (1)  can  be  expressed  in  integral  form  as: 


(23) 


(24) 


Jv-  | pH3  [g!-G2]  VpjdA  +  |  (^-^^(pHGaJdA  +  Jxt^  [pHG4]dA  = 


(25) 


Using  the  divergence  theorem,  equation  (25)  can  be  rewritten  as 


|«f-  |pH3  [o,-G2]  VpjdS  +  J  [pHG3jdA  +  ^Xi  ~  (pHG4) 


dA*0 


(26) 


Equation  (26)  integrated  over  a  control  volume  element  dA  on  the  X-Y  plane  (Figure  3a)  yields: 


pH3  [O,  -G2]  |  2  AY-pH3  [g,  -G2]  |  ,  AY  +  pH3  (g,  -G2]  JL  |  4  AX- 

pH3  |Gi -G2  ]  JL  |  j  AX+  [Xj-X,  j  JL  [pHGj  ]  AXAY  +  Xt  [pHG4]  AXAY  *0  (27) 

where  the  numeral  subscript  of  the  vertical  lines  denotes  the  intermediate  point  between  the 
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nodal  point  of  the  calculation  and  the  neighboring  points.  The  central  difference  scheme  was 
used  to  discretize  equation  (27),  however,  the  last  two  terms  in  the  equation  were  discretized 
using  the  backward  difference  scheme.  This  provided  a  stable  solution  for  the  heavily  loaded 
case.  The  energy  equation  (11)  expressed  in  integral  form  is: 


or  in  vector  form. 
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dS 

_dA 

\>- 

V 

1  1 

ij  1 2 

n  =  -eY,  0«f=-V 


VT="  §+  aj 

and  dA  =  A  H  d£  d£. 


(32) 


For  surface  5, 


n  =  -4(^e1+YCY-eZ|),  0  n  =  -  i(W-UP- Vtf 
o  o 

71  a  U  mT  9C  7  38  AH3C}  AH3C' 

and  dA  =  Vi  +^2  Vl  +  y1  <£d8. 


(33) 


For  surfaces  2,  4  and  6,  the  expressions  for  n.  If*  if  and  if  -  V  T  are  the  same  as  those  at  surfaces  1, 
3,  and  5  respectively  except  the  signs  are  reversed. 

The  volume  of  a  control  volume  element  is: 


d V  =  A  H  d£  d5  d£ 


(34) 


The  resulting  energy  equation  of  a  control  volume  in  discretized  form  is  given  by: 


-(pUT-K^)AH  |  !  A5At  +  (pUT-K-£I)AH  |  2  ASA; 


4$ 


-(pVT-K^-)AH  |  3A4AC  +  (pVT-K^-)AH  |  4  && 


a 

a 


a2  AT 


(W  -  UjS  -  Vy)pT  -  K-^j-  —■  •Vl+p2  yj  1+f  |  5  A^A5 
AH  Al 

(W-Up-VyjpT-K^j--^-  y/l+f  I  «  A§A8 

AH  At, 


=  S 


(35) 


where 
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S  =  Kp|I  |  ,  A5AC-Kg|I  |  2  ASaC  +  Ky|^  |  3  |4  ^ 

+  f  a/TT?  |  s  ^s--^  VTTpr  VTT?1"  | .  Aft  as 

+  S'  (36) 

The  left-hand  side  of  the  equation  (35)  represent  the  convective  and  diffusive  heat  transfer. 
They  are  approximated  by  the  power-law  scheme  developed  by  Patankar  (1980).  The 
discretized  form  of  equation  equation  (35)  is: 

apTp  —  %  Tg  —  awTw  —  3nTn  —  asTs  —  a^Ty  -  agTg  -  S  =  0  (37) 

where,  %  =  D2  0,  (l  -0.1  —  I  )5  +  [-F2.  ol 

aw  =D,  o,(l -0.1  |  |  )5  +  [Fi-°] 

aN  =°4  0,(1 -0.1  |  |  )5  +  [~F4.o] 

as  =D3  0,(1 -0.1  |  -^  |  )5  +  [F3,0] 

3t  =D«  0,(1 -0.1  |  ^-  |  )S  +  [-F«>°] 

aB  =DS  0.(1 -0.1  |  ^  |  )s  +  [Fs*°] 

and  3p  =  3g  +  3\v  +  3[>j  +  8s  +  3t  +  3g  • 

The  bracket  [ ,  ]  indicates  that  the  larger  one  of  the  arguments  is  selected. 
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The  expressions  for  D’s  are 


D  *** 

D'  *  45 

|  iASAC 

|  3  a$a; 

(38) 

D>-K^r 

|iiW 

The  expression  for  the  F’s  are: 

Fi  =  pUAH  | 

lASAC, 

F3  =  pVAH  | 

aAW. 

(39) 

Fj  «  jw-UfJ-Vyj^l  +  P*  Vl  +  7  |  jAJjAS. 


The  expressions  for  the  D’s  and  F’s  on  the  surfaces  2,4,  and  6  are  identical  to  those  on  the 
surfaces  1,3  and  5  respectively. 
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4.  Numerical  Procedure 


Generally,  to  obtain  a  numerical  solution  to  the  problem  of  compressible  thermal  EHD 
lubrication  of  rolling/sliding  contact,  the  discretized  Reynolds  and  the  elasticity  equations  were 
solved  using  an  iterative  technique.  The  iterative  technique  required  an  underrelaxation  factor  to 
obtain  an  accurate  and  converged  solution.  However,  for  the  heavily  loaded  conditions,  this 
underrelaxation  factor  was  usually  very  small,  which  caused  the  rate  of  convergence  of  the 
numerical  scheme  to  be  extremely  low. 


In  order  to  reduce  the  computing  time  and  increase  the  accuracy,  the  multilevel  multigrid 
technique  was  employed  in  this  analysis.  The  multigrid  procedure  used  is  similar  to  that 
implemented  by  Lubrecht  (1987)  for  the  isothermal  EHD  lubrication  of  point  contacts.  In  order 
for  the  multigrid  method  to  be  successful  it  is  essential  to  use  a  relaxation  scheme  which  has  a 
good  error  smoothing  property.  In  this  investigation,  since  the  Reynolds  equation  is  nonlinear, 
point  Gauss-Seidel-Newton  relaxation  scheme  was  used.  The  iterative  formula  for  the  pressure 
is: 


<pyrw-<pyr“ 


- ^ - (40) 

dfjj  3fij  dfkj  ,  frjj  ,  frio  dHjj  +  frij 

9Pij  '  3Pij  +  3Pij  +  dPU  ^>j-r  ®*iJ  9H‘J  ®Vi 


where,  fjj  is  the  residual  of  the  Reynolds  equation  at  node  (i,j)  and  co  is  the  underrelaxation 
factor.  The  subscript  (i+y)  denotes  the  functional  value  at  the  midpoint  between  i  and  i+1. 

Similarly,  (j+i)  also  denotes  the  functional  value  at  the  mid  point  between  j  and  j+1. 

The  energy  equation  is  also  relaxed  with  the  point  Gauss-Seidel-Newton  relaxation 
scheme.  The  viscosity  in  the  shear  heating  term  is  an  exponential  function  of  the  temperature, 
therefore,  the  change  in  energy  equation  with  respect  to  viscosity  is  included  in  the  relaxation 
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scheme.  The  iterative  formula  for  the  temperature  is: 


(T^r^cr.^r- 


^gij.k  ^&j,k 


(41) 


^Tj.k  <hlij.k 

where  gg*  is  the  residual  of  the  energy  equation  and  (a  is  the  underrelaxation  factor.  The 
complete  description  of  the  multilevel  multigrid  technique  can  be  found  in  numerous  literature 
(e.g.  Brandt,  1984).  Therefore,  the  explanation  of  the  multigrid  procedure  is  not  presented  here. 


Figure  4  illustrates  the  flowchart  for  the  computational  procedure  that  was  used  for 
calculating  the  pressure,  film  thickness  and  the  temperature  distribution  within  the  lubricant  film. 
First,  the  Hertzian  pressure  distribution  is  used  as  an  initial  guess  for  the  Reynolds  equation. 
With  the  viscosity,  density  and  film  thickness  calculated,  the  discretized  Reynolds  equation  is 
relaxed.  The  H„  value  is  modified  at  each  relaxation  to  reduce  the  residual  of  the  force  balance 
equation.  After  the  Reynolds  equation  relaxation  is  finished,  the  discretized  energy  equation  is 
relaxed.  The  ambient  temperature  is  used  as  the  initial  guess  for  the  temperature.  The  pressure, 
film  thickness  are  fixed  during  the  relaxation  of  the  energy  equation.  Then  the  pressure,  film 
thickness  and  temperature  are  transferred  to  the  next  level  according  to  the  prefixed  multigrid 
cycle  path.  The  2V  cycle  and  Full  Approximate  Scheme  was  incorporated  in  the  multigrid  path. 
The  residuals  of  the  Reynolds  equation,  the  film  thickness  equation,  the  force  balance  equation 
and  the  energy  equation  are  also  transferred.  The  film  thickness  is  much  smaller  than  the  other 
dimensions,  therefore,  semicoarsening  was  used  to  transfer  the  residuals  and  the  field  variables. 


5.  Interactive  Computer  Graphics 

An  interactive  computer  graphics  package  was  developed  for  use  in  the  design  and 
analysis  of  EHD  lubrication  of  rolling/sliding  circular  and  elliptic  contacts.  The  graphics 
package  provides  a  user  friendly  interactive  computer  graphics  interface  with  the  multigrid 
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Calculate  vucooty  ad  deaocjr 


Figure  4.  Flow  chart  for  thermal  EHD  solution  at  level  one. 
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multilevel  EHD  numerical  model.  The  graphics  package  allows  the  user  to  input  the  boundary 
conditions  and  performance  parameters  (i.e.,  load,  speed,  material  parameter,  etc.)  using 
graphics  pictorial  representation  of  the  problem.  Once  an  EHD  model  has  been  analyzed,  the 
graphics  package  enables  the  user  to  critically  examine  the  results.  It  provides  an  environment 
where  the  user  can  view  several  results  (pressure,  temperature,  etc.)  at  once.  This  allows  the 
critical  examination  of  the  effects  of  various  performance  parameters  on  the  output  results. 

Figure  5  illustrates  two  spheres  in  contact,  the  input  for  the  two  spheres  and  the  properties 
of  the  lubricant  separating  the  surfaces.  The  properties  of  the  solids  or  the  lubricant  can  simply 
be  changed  by  placing  the  arrow  with  the  mouse  cursor  on  the  appropriate  sphere,  lubricant,  etc. 
and  pressing  the  mouse  button.  Figures  6  and  7  depict  other  geometries  which  the  EHD 
lubrication  model  can  analyze.  Figure  6  illustrates  two  ellipsoid  in  contact  where  the  radius  of 
curvature  in  the  rolling  direction  (x)  is  smaller  than  in  the  y-direction.  However,  Figure  7 
demonstrates  the  condition  where  the  radius  of  curvature  in  the  rolling  direction  is  larger  than  in 
the  y-direction.  Note  that  the  ellipticity  ratio  for  the  conditions  considered  in  Figures  6  and  7  are 
1.581  and  0.6325,  respectively.  Figure  8  depicts  the  solid  properties  for  two  spheres  in  contact. 
In  order  to  change  any  of  the  material  properties,  the  cursor  is  placed  in  the  appropriate  box  and 
the  mouse  button  pressed.  A  prompt  appears  and  the  change  can  be  entered.  Every  time  any  of 
the  properties  of  the  solid,  lubricant,  etc.,  are  changed,  the  program  updates  the  geometry  and  the 
nondimensional  parameters  to  be  used  in  the  EHD  numerical  model.  When  all  of  the  changes 
have  been  entered  to  the  default  data  base,  the  EHD  model  can  be  used  to  analyze  the  input  data 
by  placing  the  cursor  in  the  "process”  menu  bar  and  clicking  on  "Run  EHD  Program."  A  box 
will  appear  where  the  output  of  the  analysis  can  be  stored.  Once  the  output  hie  has  been  entered, 
the  computer  aided  design/graphics  program  prompts  the  EHD  numerical  model  to  run  the 


Figure  5.  Geometry,  Material  and  Lubricant  Selection  for  EHD  Analysis,  Two  Spherical 
Bodies  in  Contact 
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Figure  6.  Geometry,  Material  and  Lubricant  Selection  for  EHD  Analysis,  Two  Ellipsoidal 
Bodies  in  Contact.  (Ellipticity  ratio  =  1.581) 
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Figure  7.  Geometry,  Material  and  Lubricant  Selection  for  EHD  Analysis,  Two  Ellipsoidal 
Bodies  in  Contact  (Ellipticity  ratio  =  0.6325) 
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Figure  8.  Material  Selection  and  the  List  of  Variables  That  can  be  Changed. 
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particular  data  get  developed. 


Figures  9  through  15  depict  the  results  for  the  conditions  shown  in  Figures  5  through  7. 
Figures  9  through  11  illustrate  the  pressure,  temperature,  film  thickness  contour  and  the 
volumetric  rendering  of  the  shear  stress  in  the  xz  direction.  Figure  9  allows  the  user  to  quickly 
examine  the  location  and  relative  magnitude  of  maximum  pressure  and  temperature.  The  user 
can  manipulate  and  rotate  these  figures  and  view  them  from  any  direction.  He/she  can  view 
contour  and  surface  plots  of  all  of  the  output  variable  (i.e.,  velocity,  shear  stress,  etc.).  The 
results  can  be  viewed  in  Metric  or  nondimensional.  Figure  10  illustrates  the  contour  plot  of  the 
film  thickness.  However,  the  user  can  view  the  contour  plot  of  any  of  the  output  results  (i.e., 
temperature,  pressures,  shear  stress,  etc.)  by  placing  the  mouse  button  on  the  "parameter"  in  the 
menu  bar  and  choose  the  appropriate  variable.  The  user  can  also  choose  the  number  of  color 
contours  (5  to  200)  in  the  menu  bar.  Figure  11  is  the  volumetric  rendering  of  the  shear  stress 
results.  In  this  case  the  iso-surface  results  for  a  range  chosen  by  the  user  in  the  lubricant  film  is 
presented.  This  will  allow  the  user  to  quickly  determine  the  magnitude  and  location  of  the 
parameter  of  interest  in  the  volumetric  region  that  the  lubricant  film  occupies. 

Figures  12  and  13  illustrate  the  pressure,  temperature  and  film  contour  for  the  condition 
described  in  Figure  6.  In  this  case  the  ellipticity  ratio  is  1.581.  The  pressure  exhibits  a  pressure 
spike  region  around  the  edge  and  the  back  of  the  contact  near  the  exit  region.  The  temperature 
contour  (Figure  12)  shows  that  the  temperature  reaches  its  maximum  in  the  center  of  the  contact. 
Figure  13  depicts  the  contour  plot  of  the  film  thickness.  Note  the  general  horseshoe  shape. 

Figures  14  and  15  depict  the  results  for  the  conditions  considered  in  Figure  7.  In  this  the 
ellipticity  ratio  is  0.6325  which  indicates  the  width  of  the  contact  in  the  rolling  direction  (x)  is 
larger  than  the  width  of  the  contact  in  the  y-direction.  In  this  case,  again  the  pressure  exhibits 
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Figure  9.  Pressure  and  Color  Contour  of  Temperature  in  an  EHD  Lubrication  of  Circular 
Contact 
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Figure  10.  Color  Contour  of  the  Film  Thickness  in  an  EHD  Lubrication  of  Circular  Contact 
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Figure  12.  Pressure  and  Color  Contour  of  Temperature  in  an  EHD  Lubrication  of  Elliptic 
Contact  (Ellipticity  Ratio  - 1.581) 
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Figure  13.  Color  Contour  of  the  Film  Thickness  in  an  EHD  Lubrication  of  Elliptic  Contact. 
(EUipticity  Ratio  - 1.581) 
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Figure  14.  Pressure  and  Color  Contour  of  Temperature  in  an  EHD  Lubrication  of  Elliptic 
Contact.  (Ellipticity  Ratio  -  0.6325) 
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Figure  15.  Color  Contour  of  the  Film  Thickness  in  an  EHD  Lubrication  of  Elliptic  Contact 
(Ellipticity  Ratio  -  0.6325) 


the  region  of  pressure  spike  and  the  temperature  is  maximum  at  the  center  of  the  contact 
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6.  Conclusion 


A  robust  and  fast  converging  numerical  technique  to  solve  the  simultaneous  system  of 
two-dimensional  Reynolds,  elasticity  and  three-dimensional  energy  equations  has  been 
developed.  The  numerical  technique,  the  coordinate  transformation,  as  well  as  the  iterative 
scheme,  have  been  presented.  An  interactive  computer  graphics  model  was  developed  to 
provide  a  user  friendly  environment  for  the  design  and  analysis  of  EHD  lubrication  of 
rolling/sliding  circular  and  elliptic  contacts.  The  package  is  fast  converging,  flexible  and 
provides  a  user  friendly  design  package  to  investigate  heavily  loaded  lubricated  contacts. 
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